Understanding highly excited states via parametric variations 
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Highly excited vibrational states of an isolated molecule encode the vibrational energy flow path- 
ways in the molecule. Recent studies have had spectacular success in understanding the nature of 
the excited states mainly due to the extensive studies of the classical phase space structures and 
their bifurcations. Such detailed classical-quantum correspondence studies are presently limited to 
two or quasi two dimensional systems. One of the main reasons for such a constraint has to do 
with the problem of visualization of relevant objects like surface of sections and Wigner or Husimi 
distributions associated with an eigenstate. This necessiates various alternative techniques which 
are more algebraic than geometric in nature. In this work we introduce one such method based on 
parametric variation of the eigenvalues of a Hamiltonian. It is shown that the level velocities are 
correlated with the phase space nature of the corresponding eigenstates. A semiclassical expression 
for the level velocities of a single resonance Hamiltonian is derived which provides theoretical sup- 
port for the correlation. We use the level velocities to dynamically assign the highly excited states 
of a model spectroscopic Hamiltonian in the mixed phase space regime. The effect of bifurcations 
on the level velocities is briefly discussed using a recently proposed spectroscopic Hamiltonian for 
the HCP molecule. 
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I. INTRODUCTION 



Determining the nature of the highly excited rovibra- 
tional states is important in the context of unraveling the 
intramolecular vibraticrad^eugy redistribution (IVR) 
pathways in moleculesaHauQErQ. At high energies the 
usual rigid rotor-harmonic oscillator descriptions are in- 
adequate due to the breakdown of the normal mode ap- 
proximation. The various normal modes get coupled 
leading to complicated spectra and patterns. The under- 
lying classical phase space usually has a mixed character 
i.e., regular regions interspersed with the chaotic regions. 
The mixed nature of the phase space gives rise to a va- 
riety of phenomena, classicafl and quantumH, which in 
turn lead to complex spectral patterns in terms of the 
intensities and splittings. 



Much of the recent workBOB'BrJEE3 on small poly- 
atomic molecules has clearly shown that the underlying 
classical dynamics plays a significant role in determining 
the nature of the highly excited states. In particular, 
the various classical phase space structures like peripdk. 
orbits, resonance zones and effects like bifurcationsEj'tZl 
leave their imprints on the eigenstates and hence provide 
important clues to IVR. Almost all of the work has been 
constrained to atmost two degrees of freedom systems. 
One of the main reasons for such a constraint have to do 
with the fact that visualizing the classical phase space 
structures, for that matter even the quantum states, is 
quite involved for systems with more than two degrees 
of freedom. Moreover, determining the relevant periodic 
orbitslia and higher dimensional objects like resonant 2- 
tori becomes difficult as the number of degrees of free- 
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dom increases. It is obvious that studying IVR in large 
molecules from such useful classical-quantum correspon- 
dence viewpoints would require one to come up with di- 
mensionality independent techniques. Although not the 
focus of the present work, it is important to note-, that 
novel dynamical phenomena like Arnol'd diffusion!^ and 
complex instabilityEj arise in systems with dimensional- 
ity greater than two. It remains to be seen wether such 
effects have important consequences for the classical- 
quantum correspondence studies of highly excited states 
and hence IVR. 

In principle, once the detailed phase space structure 
is at hand, the required information can be extracted 
by correlating with a phase space representation, of the 
quantum states, like the HusimiEil or the WigneiEj distri- 
butions. However, given that a detailed classical phase 
space |Stud«-| pf .| ey en three degree of freedom systems is 
sparsec3ESE3lH3E3 suggests that the straightforward ap- 
proach is far from easy. Ideally one would like a method 
that yields the required information in a basis indepen- 
dent fashion without the need to explicitly determine the 
classical periodic orbits, to visualize the surface of sec- 
tions and then comparing to the relevant phase space 
distributions. The search for such a method, if it exists, 
seems to be rather difficult at the present moment. In 
this work we attempt to provide a first step towards such 
an important yet difficult goal. 

The above mentioned factors clearly indicate that one 
needs algebraic measures that would unambigously pro- 
vide the information on the nature of the eigenstates. 
To this end certain measures, time dependent and inde- 
pendent, have been proposed in tho-literature. For in- 
stance, inverse participation ratioso (IPR) or dilution 
factors can yield information about the nature of the 
localization of the eigenstates. However IPR is a ba- 
sis dependent quantity and cannot provide information 
about the invariant classical structures relevant to the 
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eigenstate. It is possible to extract useful information 
by computing the IPR in diffecept basis which are opti- 
mal and physically motivated™^!. An explicit expression 
for the probability [distribution of the IPR has been de- 
rived very recentlycj. From a time-dependent point of 
view valuable information can be. obt ained by comput- 
ing various survival probabilitiea§E2lE3. Insights regard- 
ing the nature of the underlying classical phase space 
have emerged from establishing the existence of corre- 
lated intermediate timescale dynamics for rovibrational 
HamiltoniansOEy. The role played by the various classical 
structures in determining the natuxe of the survival prob- 
abilities is well knowrS. Recently^ it has been proposed 
that a dip in the level spacings is a characteristic spectral 
pattern indicative of a seperatrix in the underlying clas- 
sical phase space. This is an important observation for 
fitting the experimental spectrum with effective Hamil- 
tonians since seperatrices are indicative of qualitatively 
different dynamical behaviours coexisting in the system. 

On the other hand, when the system Hamiltonian 
is known, a number of authors have focused atten- 
tion on parametric variation of the eigenvalues of the 
systemHESEH. In particular, plots of the energy eigen- 
values versus a perturbation parameter show very spe- 
cific patterps for resonant systemsLjO. It is also well 
knownEaSEll that the parametric variation of the eigen- 
values of a system can be mapped onto a classical Hamil- 
tonian system with the eigenvalues corresponding to the 
positions of a set of fictitious particles and the parame- 
ter playing the role of a psuedotime. Within this map- 
ping the slopes of the energy curves correspond to level 
velocities!^. Thus one associates a level velocity with 
every eigenstate of the system. Significant progress has 
been made in|-studying the statistical— properties of the 
level velocitiescil and the accelarationsc3 (curvatures). In 
analogy with the random matrix theory (RMT) results 
for eigenvalue statistics, the level velocity and [Curvatures 
exhibit universal behaviours in the RMT limito. For ex- 
ample, in RMT the level velocities are found to be gaus- 
sian distributed and a relational has been established be- 
tween the IPR distribution and the level velocity distribu- 
tion. Accordingly, deviations from the RMT limit result 
in nonunivosaal distributions signifying localization of the 
eigenstatesac3'C3'El The preceeding discussion suggests 
that it might be useful to explore the possibility of us- 
ing the level velocities to understand the highly excited 
states of an isolated molecule. 

The purpose of this paper is to show that the level 
velocity of an eigenstate is strongly correlated with the 
phase space nature of the eigenstate. It will be demon- 
strated that the magnitude and relative signs of the level 
velocities give a direct indication of the important clas- 
sical structures that influence the corresponding eigen- 
states. Theoretical arguments for the correlation and nu- 
merical support are provided in the next section for clas- 
sically integrable, single resonance systems. In section 
III the level velocity "spectrum" of a nonintegrable, mul- 
tiresonant system is studied. Combined with an earlier 



work on the fingerprints of resonances on the eigenlevel 
dynamics it is possible to have an approximate dynam- 
ical assignment of the eigenstates when the underlying 
phase space is mixed. We conclude with a discussion of 
the major advantages and disadvantages of our approach 
and briefly discuss the effects of bifurcation on the level 
velocities. 



II. THEORETICAL BACKGROUND 

In this section we introduce the level velocity spectrum 
and provide arguments for the observed correlation be- 
tween level velocities and the phase space nature of the 
corresponding eigenstates. The arguments presented in 
this section are strictly applicable only to single resonant 
Hamiltonians. Nevertheless, we expect the arguments to 
be valid for the multiresonant systems exhibiting mixed 
phase space dynamics. 

Consider the following n : m resonant Hamiltonian: 



H(r) = £ ("i N i + 

.3=1,2 



+ T 



>.,/-V; 



h.c 



(1) 



with Nj — (rij + 1/2) and (aj,aj) denoting the standard 
destruction and creation operators for the j th mode. The 
eigenvalues and eigenstates of the Hamiltonian are de- 
noted by E(t) and \a(r)) respectively. The form of the 
diagonal part of H(t) is not restricted to the choice made 
above and can include higher order anharmonicities as 
well. The form of the Hamiltonian is^that of an effec- 
tive or spectroscopic Hamiltonian and Hq corresponds la 
the Dunham expansion. Defining the polyad operatorca 
P = (n/m)ni + n,2 it is easy to see that [H , P] = 0. 
Thus H is block diagonal with each block labeled by the 
polyad number P and all the eigenstates can be assigned 
in terms of two quantum numbers as |a(r)) = \P, v(t)). 
The cumiitum number v — 0,1,... is an excitation 
indexo'EZI whose range is related to the strength of the 
resonant coupling r. The eigenstates for a given P are 
linear combinations of the zeroth order basis functions 
n) = \ni,ri2) with (njra)ri\ -\-n% = P. Hence we can 

write |a(r)) = SneP c n"( T )l n ) with c n»( T ) = i n \ a ( T ))- 
The level velocity associated with an eigenstate |a(r)) 
is obtained using the Hellman-Feynman theorem and 
given by: 

x a (r;P) = = (a(r)\V n:m \a(r)) (2) 

with V n:m being the perturbation. It is possible to write 
the level velocity in terms of the expansion coefficients 
Cnaij) as 



n.n'eP 

c q (t;P)= c* na ( T )cn> a (T)V n 



(3) 
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where the matrix element V nn i = (n|V^j :m |n'). At this 
stage we summarize the salient features of x a (r;P) for 
single resonant systems. 

1. The average of the level velocities for states belong- 
ing to a given polyad P vanishes irrespective of the 
coupling strength r. This can be easily seen since 

n.n'eP 

J2i a (T;P) = fn»'Z>'l J »( J >l n > 

n.n'eP 

= J2 V m ,(n'\n) = (4) 



This fact gives rise to the "fa: 
the energy correlation diagra™^ 
tion parameter being r. 



" structures in 
with the varia- 



2. The leszpl velocities exhibit linear parametric 
motiorO beyond a certain critical coupling 
strength r c . The critical coupling-can be approxi- 
mately determined via a-ChirikovO analysis. It has 
been argued previously^ that the asymptotic ve- 
locities for a 1 : 1 state |P; v) are given by [P — 2v). 
In general, one does not expect the resonant cou- 
pling strength in a molecule, to be close to r c . Con- 
sequently the velocities, except for states localized 
at or near the center of the resonance zone, do show 
nonlinear parametric variation. However, as an ex- 
ample we mention that the fermi resonance cou- 
pling strength in the spectroscopic effective Hamil- 
tonian of the CHF-p molecule does come close to 
the critical valueca. For a generaJ_jn : m reso- 
nant Hamiltonian it was conjecturedEZI, and numer- 
ically confirmed, that the velocity of the state \P; 0} 
scales with the polyad as p(« l + n )/ 2 . 

It is clear that through the velocity one is studying the 
response of the eigenstate to a specific perturbation. It is 
not surprising that states localized in or near the specific 
resonance zone would have the largest response. That 
the response is strongly correlated to the phase space na- 
ture of the eigenstate is interesting and we now provide 
plausible reasons for such an observation. We note that 
there has been an earlier work wherein similar correla- 
tions were observedEJ. 

To start with note that the matrix element 



(n|(4) m (a 2 )"|n') = V nm (n[ , n' 2 )S ni tn > i+m 8. 



?i2 -n ~n 



(•5) 



yields an expression for the level velocity in the following 
convenient form: 



p 



x a (r; P) = 2 ^ K m «; PJ^^^t) (6) 

n{=0 

where we have supressed the second index on the coef- 
ficients. From a classical-quantum correspondence per- 
spective it would be ideal to have a semiclassical expres- 
sion for x a (r;P). The shifted overlap like structure in 



the above expression for-J; ct (r; P) is strongly reminiscient 
of the Wigner functioned. To make the correspondence 
explicit the level velocity is written as: 



x a (r;P) 



(a\V n:m \a) 
Tr 



(7a) 
(7b) 



V n:m \a}{a\ 

dId<f>V w {I,4>)W a {I,4>) (7c) 



1 

~ 2vrft 

where W a (I,4>) is the Wigner function associatedLjHath 
the state |a) and Vw{I,<j>) is the Weyl symboOtBEj of 
the operator V n:m defined by 



b > e -W/K^ + <l//2\V nsm \<t> - <t>'/2) (8) 



Note that the Wigner function is the Weyl symbol of 
the density operator. Although the system is two dimen- 
sional it is sufficient to focus on an effective one dimen- 
sional Hamiltonian due to the existence of the polyad. 
Thus (/, cj>) are action-angle variables pertinent to the 
classical limit of the reduced zeroth order Hamiltonian 
(cf. appendix). In particular, all of the expressions in 
this section are for a given polyad P. This already sug- 
gests the role pWed by classical phase space structures 
since it is knownO that for integrable systems the Wigner 
function condenses onto the invariant torus correspond- 
ing to the eigenstate. 

In order to determine the Weyl symbol Vw (I, 0) con- 
sider 



if e -^'/ n (0 + ( /)72|(a t 1 ) m (a2) n |0-072) 



E- 

^ 2?r 



i —ilcji' /h i(m— ni)</>gi(ni+ni)<£'/ 2 



X)Wni;P). 



,im4> 



(9) 



Sinc^j-the Weyl symbol of a hermitian operator is 
reaSo, Vw{I, 4>) is given by: 

V W (J,4>) =2V nm (I/fi-m/2;P)cosm(f> (10) 

The Wigner function is now approximated!^] in order to 
elucidate the role of the fixed points. First the wavefunc- 
tion is written in its primitive semiclassical form 



Mi) ~J2 Ab *W e ~ 



iS ba {i)/h 



(11) 



where the index b denotes the branch. The phase and 
amplitude are given by 



S ba (I) 
Aba 



/I ~ 
fba(l')dl' + Y^b a (12a) 

dH cl (I 



-1/2 



(12b) 
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The Maslov index0 is denoted by fiha an d 
H c i(ifiba(I) , I) — E a . Next the Wigner function is 
approximated as 



W a (I,<f>) = /dcre to */Va(/ + ff/2)^*(/-ff/2) 

« 27r^|A &ct (/)| 2 <5[</>-< / ? bQ (/)] (13) 



by expanding the amplitude and phase terms to first or- 
der. The intrabranch terms have been neglected in the 
approximation. The above result for the Wigner func- 
tion highlights it as the classical density on an invariant 
torus in the corresponding phase space. Combining the 
approximate Wigner function with the Weyl symbol for 
Vn-.m one arrives at the central result of this section- a ap- 
proximate semiclassical expression for the level velocity 
associated with an eigenstate \a) as 

f Pc 1 

i Q (r;P)«2V/ dI-^V(I;P c )cos[m^ ba (I)} (14) 
b Jo \I\ 

where P c is the classical analog of P and the leading order 
Ji term of V nm {I /h-m/2; P) is denoted by V(J; P c ). The 
above semiclassical expression emphasizes the fact that 
fixed points of the classical Hamiltonian i.e., 1 = play 
an important role in the level velocities. Moreover for a 
single resonant Hamiltonian the phase term in the level 
velocity expression is ±1 depending on the nature of the 
fixed point. Thus, relative signs of the level velocities 
are important. In the appendix further insight is pro- 
vided on the polyad scaling of the velocities at such fixed 
points by a classical analysis of the function V(I;P C ). 
Based on the above analysis it is expected that the level 
velocities will be large in magnitude for states localized 
around the various fixed points. A_similar expression for 
the level velocity can be derivedE3 independently start- 
ing from the Eq. (^j) and following an approach based on 
psuedodifferential equations^. 

To provide support for the theoretical arguments pre- 
sented above consider an effective 2-mode Hamiltonian 
H = Hq + tVi-x which is 1 : 1 resonant with the pa- 
rameters u)\ = 1,W2 = 0.8, an = —0.03, a22 = —0.02, 
and r = 0.01. In order to determine the phase space na- 
ture of the various eigenstates the Husimi distributions 
are computed. Although one could compute the Wigner 
functions associated with the states the Husimi functions 
are preferred due to their smooth natural. 

In Fig. ([!]) we show the level velocities for the states 
belonging to polyad P = 8. The three states which lie 
inside the 1 : 1 resonance zone are indicated by the label 
v. It is easy to see that the coupling is below the critical 
threshold since the velocity of the 8; 0) state is less than 
the asymptotic value of 8. In Fig. (||) the Husimis as- 
sociated with the states are shown superimposed on the 
appropriate Poincare surface of sections. For the sake of 
clarity only the contours with the maximum values are 
shown. The state 8; 0) is localized on the stable periodic 
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FIG. 1: Level velocity spectrum associated with the states 
belonging to the polyad P — 8 for the single 1 : 1 resonant 
Hamiltonian. 



orbit i.e., at the center of the 1 : 1 island. As one pro- 
ceeds to the states with higher v the Husimi maximum 
moves out of the island center and for the state |8;2) 
eventually localizes on the unstable periodic orbit i.e., 
around the seperatrix. It is important to note that the 
state 1 8; 0) has a large positive velocity and the state |8; 2) 
has a large negative velocity. This observation is consis- 
tent with the results of the theoretical analysis. The rest 
of the states in P = 8 are either local modes or distorted 
local modes as evident from the Husimis in Fig. (0). 

In order to further check the correlation between level 
velocities and the phase space nature of the correspond- 
ing states the velocities and Husimis are computed for 
an integrable 2 : 1 system H = Hq + rV^a with the 
coupling strength r = 0.01. The level velocities for the 
states belonging to the 2 : 1 polyad P = 14 are shown in 
Fig. (|^). Again it is clear that the coupling is below the 
critical coupling and the level velocities show a charac- 
teristic pattern. The state 1 14; 0) is expected to localize 
around the 2 : 1 stable periodic orbit whereas the state 
1 14; 3) should be localized around the unstable periodic 
orbit. The husimis shown in Fig. (^) essentially confirm 
the expectations. In both the integrable cases, for cou- 
pling strengths near or above the critical value, the level 
velocities are large positive for the state in the center of 
the resonance island and equally large but negative for 
the seperatrix state. Thus the relative signs of the level 
velocities are crucial and provide information on the fixed 
points that dictate the phase space nature of the eigen- 
states. The absolute sign is clearly unimportant since 
a change in the sign of the coupling would interchange 
the signs of the velocities. In addition the state localized 
about the stable periodic orbit exhibits linear parametric 
motion and the velocity is closer to the classical value, 
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FIG. 2: Husimi distributions associated with the states be- 
longing to the polyad P — 8 for the single 1 : 1 resonant 
Hamiltonian. For comparison the Husimis are superimposed 
on the surface of sections. The top three panels show the 
resonant states with v — 0, 1, 2. 



FIG. 4: Husimi distributions associated with the states be- 
longing to the polyad P = 14 for the single 2 : 1 resonant 
Hamiltonian. The top three and the middle left panels corre- 
spond to the resonant states with v = 0, 1, 2, 3 respectively. 
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FIG. 3: Level velocity spectrum associated with the states 
belonging to the polyad P = 14 for the single 2 : 1 resonant 
Hamiltonian. 



computed using Eq. ( 27b ), as compared to the state lo- 
calized about the unstable periodic orbit. 

The results of this section clearly indicate that the level 
velocities associated with the eigenstates are correlated to 
the structures in the classical phase space. In the follow- 
ing sections additional numerical examples are provided 
which highlight the correlation and the potential for us- 
ing the level velocities to assign highly excited states of 



multiresonant Hamiltonians. 



III. APPLICATION TO NONINTEGRABLE 
SYSTEMS 

The states of single resonant Hamiltonians studied in 
the previous section can be assigned in various 
Complications due to bifurcations can be sorted outU^ 
with a detailed study of the classical-quantum correspon- 
dence for the system. However with the addition of an- 
other independent resonance the polyad constant of the 
motion is destroyed and assignments of the states be- 
comes a nontrivial task. The underlying classical phase 
space can exhibit a rich variety of dynamics ranging from 
the near-integrable to the mixed to full chaos. An im- 
port anL-ftJasGijatkua- that has emerged from the various 
studiesli300'oL§li3 is that classical phase space struc- 
tures can form the basis for a dynamical assignment of 
the various eigenstates of the system. The advantage of 
such a dynamical assignment is that it is based on invari- 
ant structures and directly reflects the relevant molecular 
motions in some energy range of interest. More impor- 
tantly, such dynamical assignments are possible over the 
entire range of the phase space dynamics. For instance, 
a completely chaotic phase space at the outset seems 
to preclude any meaningful assignment of the states. 
Even for such hard chaotic systems it is wellknown that 
the eigenstates can.hij localized due to various quantum 
and classical effect soo'ElotH A prime example is the 
scarringL3 of the eigenstates by .the periodic orbits which 
has been studied in great detailEa over the years. In this 
case the relevant periodic orbits provide an avenue to 
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dynamically assign the eigenstates. 

From the perspective of studying IVR in isolated 
molecules a key observation is that the underlying phase 
space overwhelmingly belongs to the class of mixed dy- 
namics. The fact that the effective spectroscopic Hamil- 
tonians used to gain insights into IVR have the struc- 
ture of a hierarchical local random matrixta is largely re- 
sponsible for the mixed nature of the phase space. Dy- 
namical assignment of the highly excited eigenstates in 
such mixed phase space regimes is important and highly 
relevant. A significant bottleneck to such a classical- 
quantum correspondence based approach arises for sys- 
tems with dimensionality greater than two. Any dynam- 
ical assignment of the eigenstates necessarily involves 
finding the important classical structures in the phase 
space and correlating them with the eigenstates or some 
phase space representation of the eigenstates. Already 
for a three degree of freedom system the Poincare sur- 
face of sectioned is four dimensional and so is the com- 
puted Wigner or Husimi distribution of the eigenstates. 
This leads to difficulties in visually correlating the states 
with appropriate classical objects. One way to overcome 
this problem is to have algebraic measures and in this re- 
gard the results of the previous section suggest the level 
velocities as one possible candidate. 

In order to assess the utility of the level velocities in 
providing a dynamical assignment of the eigenstates a 
model two resonance effective Hamiltonian H = Hq + 
TliVi-i + 721^2:1 is chosen. The parameters of the ze- 
roth order Hamiltonian are as given in the previous sec- 
tion. The resonant coupling strengths are chosen to be 
Tn = T21 = 0.01. The Hamiltonian is nonintegrable and 
the energy range [8, 9] is chosen as an example since the 
phase space exhibits mixed dynamics. The level veloc- 
ities x Q (rn) and x a (r2i) are computed for some of the 
eigenstates in the selected energy range. Note that in the 
nonintegrable case the velocities can be expressed as: 

p 

(rjjlS) 

p m=o 

= 2^x Q (r;P) (16) 
p 

In this work the x a (r;P) will be called as the partial 
level velocities. The level velocities are scaled as 

x q (t) = — (x q (t) - (x q (t)}) (17) 

OE 

where oe is the variance of the velocities and (x a ) is the 
average of the level velocities. This is done since the ve- 
locities scale differently with the polyad P and we would 
like to compare the two velocities for a given eigenstate. 
The scaled velocities are dimensionless, zero centered and 
have unit variance. 

In Fig. ([5]) we show the scaled level velocity spectrum 
for selected eigenstates in the energy range [8, 9]. In Ta- 
ble I the state energies, state IPRs in various basis, the 



TABLE I: Assignment of selected states for the nonintegrable 
system. 



State E Lq Lu L21 in ±21 Assignment™ 

68 8.038 0.35 0.39 0.89 -0.17 1.08 P21 - 1 

69 8.044 0.84 0.85 0.98 0.01 -0.79 dhi 

70 8.164 0.22 0.28 0.79 -0.30 -1.56 u,P 21 

71 8.189 0.74 0.95 0.77 -0.65 -0.08 din 

72 8.203 0.99 1.00 0.99 -0.03 0.13 I 

74 8.253 0.25 0.28 0.71 -0.53 -0.08 * 

75 8.277 0.35 0.37 0.85 -0.05 2.58 s,P 2 i-l 
77 8.313 0.21 0.75 0.26 -1.54 -0.38 u,Pn 

79 8.338 0.94 0.99 0.95 -0.15 0.03 I 

80 8.367 0.33 0.70 0.41 1.64 -0.46 P n 
82 8.423 0.99 1.00 0.99 -0.02 0.13 I 

85 8.460 0.36 0.71 0.43 3.81 -0.38 s,Pu 

86 8.473 0.21 0.23 0.81 0.11 1.07 P 2 i 

88 8.514 0.38 0.39 0.93 0.08 -1.34 u,Pn+l 

89 8.583 0.99 1.00 0.99 -0.01 0.13 I 

91 8.650 0.83 0.97 0.85 -0.41 -0.08 din 

92 8.678 0.95 0.99 0.96 -0.11 0.03 I 

93 8.682 0.22 0.25 0.76 -0.40 -0.27 * 

94 8.683 0.99 1.00 0.99 -0.01 0.13 I 
96 8.719 0.33 0.35 0.85 -0.05 2.89 s,P 2 i 

98 8.763 0.48 0.47 0.96 -0.03 -1.33 u,P 2 i+2 

99 8.848 0.24 0.58 0.27 -2.24 -0.38 u,Pn+l 

100 8.869 0.18 0.20 0.27 -0.39 -0.17 * 

101 8.898 0.18 0.19 0.53 0.35 0.84 *,2 : 1 

102 8.958 0.96 0.99 0.97 -0.09 0.03 I 

103 8.967 0.15 0.39 0.21 0.93 -0.83 *,l:l,Pii + l 

104 8.972 0.37 0.40 0.51 0.24 -0.92 *,2 : 1 



a s, u, I, and dl denote stable, unstable, local, and distorted local 
respectively. States involved in avoided crossings are marked with 
a *. See the text for details. 

scaled level velocities and the assignments are provided. 
The IPR of a state \a) in a basis \b) can be obtained as 

L,(a)=]T|(&|a}| 4 (18) 

b 

In Table]!], Lq, Lu, and L21 denote the IPRs in the zeroth 
order number basis, the 1 : 1 basis and the 2 : 1 basis 
respectively. The standard interpretation of the IPR is 
that a low value signifies extensive mixing in the basis of 
choice whereas a high value indicates very little mixing of 
the basis states. The scaled velocity spectrum in Fig. (||) 
clearly shows different classes of states for which a brief 
description is provided below. 

A. Local mode states 

These states (72,79,82,89,92,94,102) are neither influ- 
enced by the 1 : 1 perturbation nor the 2 : 1 perturbation 
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State # 



FIG. 5: Scaled level velocity spectrum associated with the 
states in the energy range [8, 9] for the multiresonant Hamil- 
tonian. The solid lines with filled circles correspond to the 
1 : 1 case and the dashed lines with open squares correspond 
to the 2 : 1 case. 



FIG. 6: Husimis for selected states superimposed on the cor- 
responding surface of sections. See the text for details. 



C. Resonant states 



and hence nonresonant. The scaled velocities i a (Tn) and 
£a(T~2i) are both very small indicating the local nature 
of these states. Note also that the IPRs are close to one 
in every basis for these set of states. As an example, the 
husimi of state 94 superimposed on the surface of section 
is shown in Fig. (|^a) . The zeroth order quantum numbers 
(ni,ri2) are good quantum numbers for the local mode 
states. In Table |such states are denoted by T. 



B. Distorted local mode states 

As suggested by the name the distorted local mode 
states (69,71,91) are states that are lying close to either 
one of the resonance zones. Thus the husimi distribution 
of such states are distorted but topologically nonreso- 
nant. One such example is shown in Fig. (Ob) for the 
state 71. The level velocities clearly indicate the partic- 
ular resonance zone to which the states are in proximity. 
This is precisely the reason that the state 71 is labeled as 
' cZZ 1 1 ' . The IPR values are also high and essentially agree 
with the assignments from the level velocity spectrum. 
The zeroth order quantum numbers are not so good a la- 
bel for these distorted states. It is also possible to think 
of these states as a kind of 'borderline' states and hence 
a slight increase in the appropriate resonance coupling 
strength will render them resonant. 



The resonant states (68,70,75,77,80,85,86,88,96,98,99) 
show up in the velocity spectrum in Fig. (|^) as states 
with large positive or negative velocities. In accordance 
with our discussions in the previous sections these states 
are localized in one or the other resonance zone. More- 
over states with large positive velocity are expected to 
be localized about the corresponding stable periodic or- 
bit. On the other hand states with large negative veloci- 
ties are expected to be localized about the corresponding 
seperatrix region. For example state 85 has a large pos- 
itive £ q (tii) and comparitively negligible i a (T~2i)- Con- 
sequently this state is anticipated to be localized exclu- 
sively in the 1 : 1 resonance zone and the Husimi distri- 
bution to be localized about the 1 : 1 resonance island in 
the underlying phase space. That this is indeed the case 
is confirmed in Fig. (J^e). Thus state 85 is labeled by V 
indicating localization about the stable periodic orbit. In 
Table | we have also appended a label P\\ suggesting that 
the particular 1 : 1 polyad Pn is almost a good quantum 
number to assign the state. Similarly state 96 is local- 
ized about the 2 : 1 island according to the level velocities 
and Fig. (^|g) supports the assignment. An approximate 
polyad P21 is also attached to this state. 

In order to complete the assignment the approximate 
polyads Pn and P21 need to be determined. To this end it 
is useful to examine the partial level velocities i a ij; P) as 
a function of P for the corresponding states. P is chosen 
appropriately for the particular resonance. Thus, P = 
ni+n2 for the 1 : 1 resonance and P = 2ni+ri2 for the 2 : 
1 resonance. It is important to note that this information 
is obtained enroute to computing the level velocity of a 
state and does not require a seperate computation. In 



8 



0' 
-0.1 
-0.2 
-0.3 
-0.4 
-0.5j 









- 


a 




■ 










4 _ 

, i , 



>- ' 

!■- 
t _ 


' e - 

l 


i _ m / 

- , i , ik, i , 




FIG. 7: Unsealed partial level velocities for selected states as 
a function of the polyad. The labeling of the figures corre- 
sponds to that in Fig. (^|). The solid lines with filled circles 
correspond to the 1 : I case (P = 711+112) and the dashed lines 
with open squares correspond to the 2 : 1 case (P — 2ni+n2). 



Fig. (ffkg) the partial velocities are shown for states 85 
and 96 respectively. It is clear from the figure that Pn = 
10 and P21 = 14. A similar analysis can be done on all 
the states and as an example the partial velocities are 
shown in Fig. (|?|a) for the state 94. One can immediately 
assign this local mode state as (m, n 2 ) = (0, 15) 

Once the resonant states with large positive velocities 
have been assigned the rest of the resonant states can 
be easily assigned as shown in Table Q. Two important 
properties of the level velocities guide the assignments. 
The first property is that for positive coupling constants 
the states localized about the stable periodic orbits with 
some approximate polyad P are energetically higher than 
the states localized about the unstable periodic orbits be- 
longing to P. Exceptions can occur due to bifurcations 
and this is discussed in the following section. Secondly, 
velocity of states with some approximate P scale differ- 
ently with P for the 1 : 1 and the 2 : 1 resonances. Thus 
states belonging to the 2 : 1 resonance are seperated more 
as compared to the, states belonging to the 1 : 1 reso- 
nance. A ChirikovE.3 like estimate is also an indicator of 
the number of states bound by the resonance. Moreover, 
studying the partial velocities helps in confirming the as- 
signments. In Figs. (||c,f) we show the Husimi distribu- 
tions for the states 77 and 70 respectively. The Husimis 
are localized in the chaotic regions corresponding to the 
integrable limit seperatricies and lend support to the as- 
signments in Table | where the label V is used to denote 
such states. 
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FIG. 8: (a) Sharp avoided crossing between states 74 and 75 
as a function of the 2 : 1 coupling parameter T21. (b) Broad 
avoided crossing between states 103 and 104 as a function of 
the 1 : 1 coupling parameter Tii. States are analysed in this 
work for m = T21 = 0.01. 



D. Mixed states 

There are a few states (74,93,100,103) which cannot 
be assigned in a straightforward fashion. As an example 
state 74 according to the IPR values should be predom- 
inantly a 2 : 1 resonant state. However the level veloc- 
ities show the opposite trend. In Fig. (^h) we show the 
Husimi for state 74 superimposed on the surface of sec- 
tion. It is clear that the state is delocalized over both the 
chaotic as well as the regular regions of the phase space 
as compared to the other states in Fig. The Husimi, 
however, is peaked in the 2 : 1 regions of the phase space 
in agreement with the IPR values. These observations 
are symptomatic of avoided crossings since at an avoided 
crossing regular and chaotic states can mix. Indeed on 
examining the variation of the eigenvalues with the 2 : 1 
coupling parameter T21, shown in Fig. (^a), it is seen 
that the state 74 is involved in an avoided crossing with 
the regular state 75. The important observation is that 
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the state 74 is located close to the center of the avoided 
crossing. This explains the low value for the velocity 
x a {T2i) and the apparent difficulty in the assignment of 
the state 74. The partial velocities for the state 74 shown 
in Fig. clearly demonstrate the delocalized nature of 
the state. The assignment of the state 93 is complicated 
due to similar reasons. The states 100 and 101 are also 
close to being involved in an avoided crossing. In this 
case the velocities can be used to provide a nominal as- 
signment. 

On the other hand states 103 and 104 are again close 
to being involved in an avoided crossing as shown in 
Fig. (bp). Note that the parameter being varied in 
Fig. (Bp) is the 1 : 1 coupling strength t\\. The husimi 
distribution for the state 103 is shown in Fig. (^i) and 
looks topologically very much like the case in Fig. (||d) 
except for the peak in the 2 : 1 region. Again the level 
velocities can be used to nominally assign the state 103 
as a 1 : 1 state belonging to the polyad P\\ + 1. Thus the 
mixed states arising in this energy range are primarily 
due to two state avoided crossings. 
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IV. DISCUSSION AND SUMMARY 

The main focus of the present work is to come up with 
techniques which would allow a dynamical assignment of 
highly excited states of multiresonant, multidimensional 
(> 2) Hamiltonians. It is well known that an intimate 
knowledge of the various classical phase space structures 
is the first important step towards such dynamical as- 
signments. However obtaining a detailed picture, akin to 
the two dimensional systems, of the phase space struc- 
ture for systems with three or larger degrees of freedom 
is diffcult. In this paper we have argued that the phase 
space nature of an eigenstate is encoded in the manner 
that it responds to variation of specific parameters of 
the Hamiltonian. In other words, the level velocity as- 
sociated with an eigenstate is strongly correlated to the 
phase space nature of the eigenstate. Theoretical argu- 
ments for the observed correlation have been provided for 
single resonant systems. In a way the current work be- 
gins to provide a basis for understanding and utilizing the 
usual energy correlation diagrams. We have shown the 
utility of the level velocities in a dynamical assignment 
of the highly excited states of a 2-mode multiresonant 
Hamiltonian. It is important to point out here that our 
approach does not require explicit determination of the 
periodic orbits for the system. It is equally relevant to 
mention that the level velocity approach is independent 
of the dimensionality of the system and is manifestly ba- 
sis invariant. Thus the current work is a first step in 
coming up with algebraic measures facilitating the study 
of the classical-quantum correspondence for multidimen- 
sional systems. 

The method proposed in this paper for the dynamical 
assignment of highly excited states, however, has a few 
weakpoints. Firstly, it is quite possible for the level ve- 



FIG. 9: Level velocity spectra for the HCP molecule for the 
various CP stretch-bend polyads P. [r] and [B] are stable 
periodic orbits, (a) P — 14 (b) P = 18. Notice the pertur- 
bation of the positive part of the level velocities. / labels the 
isomerization state, (c) P = 22 (d) P = 32. 

locity of a particular eigenstate to be close to zero. For 
example, for coupling strengths near the critical value, 
the 1 : 1 state |6;3) will exhibit zero velocity. Assign- 
ing such states solely on the basis of level velocities will 
be problematic. Secondly, as demonstrated in the last 
section, states involved in an avoided crossing can pose 
problems for the assignment as well. This is not surpris- 
ing since such states arc difficult to assign by any of the 
existing methods. Nevertheless, taking linear combina- 
tions of the two states and assigning the demixed states 
is one possible solution. In both instances it is still pos- 
sible to gain insights into the states by examining other 
measures like the IPRs in conjunction with the level ve- 
locities. 

One of the important issues that we have not addressed 
has to do with the effect of the various bifurcationsli3ll3'EZl 
on the level velocity spectrum. In this work we have 
explicitly shown the important role played by the fixed 
points in the classical phase space as far.as the level ve- 
locities are concerned. BifurcationsESEa are associated 
with creation or destruction of fixed points and hence 
signal qualitative change in the dynamics. Thus, it is 
clear that the level velocities would be quite sensitive 
to bifurcations. Work is in progress in our group to 
study the fingerprints of bifurcations on the level ve- 
locity spectrum in a sytematic fashion. However, we 
conclude this paper by showing some preliminary work 
on the HCP molecule in order to gauge the utility of 
the level velocities. We have chosen HCP as an exam- 



10 



pie due to the fact rthat the system has been investi- 
gated experimentallyEd and thoroughly studied from the 
classical-quantum persnectiveMEJ. The effective spec- 
troscopic HamiltonianEJ for HCP involves a 1 : 2 Fermi 
resonant coupling between the bend and the CP stretch. 
Thus there exists a polyadEZI P = V2 + 2v% with v 2 and 
«3 refering to the bend and the CP stretch modes re- 
spectively. The states can be assigned in terms of the 
vibrational angular momentum, the number of quanta in 
the CH stretch («i), and the polyad P. Further for ev- 
ery given P there are P/2 + 1 levels indexed by i = 
at the top of the polyad and i = P/2, at the bottom of 
the polyad. It has been cstablishcdL3 that for v± = 
and up to P — 14 the states evolve smoothly from be- 
ing pure CP stretch at the bottom of the polyads to be- 
ing approximately pure bend at the top of the polyads. 
However starting at P — 18 a new family of states, 
called isomerization states, come into existence due to 
a saddle node bifurcation of the periodic orbits. Above 
P = 22 the isomerization states are almost pure bend 
and closely follow the minimum energy path leading from 
HCP to CPHEj. Keeping this information in mind in 
Fig. (j^) we show the unsealed level velocity spectrum cor- 
responding to polyads P = 14, 18, 22, and 32 for the HCP 
molecule. It is crucial to note that for the effective spec- 
troscopicJIamiltonian fit the Fermi resonant coupling is 
negativeE3. We thus anticipate a large negative velocity 
for the state localized around the stable periodic orbit. It 
is also known that for low polyad values there are two sta- 
ble periodic orbits denoted by [r] and [£>]. The velocity 
spectrum for P = 14 shown in Fig. (||a) clearly indicates 
that the i = state is localized around the [B] periodic 
orbit. From our theoretical analysis and the appendix we 
expect the state associated with the [r] periodic orbit to 
have a very small velocity which is confirmed in Fig. (^a) 
for state i = 7. The velocity spectrum suggests that most 
of the states are resonant and indeed the classical esti- 



mate (-28.51) (cf. Eq. fl27bD ) for the velocity of the state 
i = is very close to the computed value (—26.36). In 
Fig. (^|b) the velocity spectrum is shown for P = 18. For 
this case we clearly see a pertubation in the spectrum in 
the positive velocity region and the resulting isomeriza- 
tion state is labeled as /. The negative velocity region 
is unperturbed. The perturbation is quite prominent in 
Figs. (^|c,d) corresponding to the P — 22 and the P = 32 
cases. In every case the velocity of the [B] state is close 
to the classical estimate. Given that the negative portion 
of the velocity spectrum is unperturbed we conclude that 
the bifurcation is creating a stable periodic orbit (SN) 
and an unstable periodic orbit (SN) at ip = ir/2 in the 
reduced phase space. Thus the isomerization states form 
a sequence from being localized about the SN periodic 
orbit to the SN periodic orbit. In this instance as the 
saddle node bifurcation creates the stable/unstable pair 
at the same angle rp = ir/2 the level velocities for states 
localized about both the periodic orbits are expected to 
be large positive. The caricature of the phasG-spsuce just 
described agrees well with the earlier studicsli§E3. It is 



clear that the level velocity spectrum is correlated with 
the phase space structure and capable of detecting non- 
trivial changes in the phase space due to the bifurcations. 

This work is a first step towards developing tools 
suitable for studying classical-quantum correspondence 
in multidimensional systems. The systems studied in 
this paper are certainly encouraging and warrant fur- 
ther study of the level velocities and different measures 
based on the level velocities-JOne such sensitive measure 
has been proposed recentljo which correlates the level 
velocities with the eigenstate resolved spectrum of the 
molecule in order to gain insights into the IVR process. 
Currently we are exploring the utility of the level veloc- 
ities in assigning the states of three degree of freedom 
multiresonant Hamiltonians which are effectively two di- 
mensional due to the existence of a superpolyad num- 
ber. Spectroscopic Hamiltonians have been determined 
for many such effectively two pdimensional systems like 
H 2 OE3, CHBrClFE, and DC(M Preliminary studiesEl 
on such systems indicate that the level velocity spectrum 
is well suited for the purpose of dynamical assignment of 
the eigenstates. Studies are also underway of true three 
degree of freedom systems in order to test the approach 
and for insights into the highly excited states of such 
systems from the classical- quantum correspondence per- 
spectives. 
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Appendix: Classical insight into the level velocities 

In this appendix we provide a classical argument for 
the scaling of the level velocities with the polyad P for a 
general single resonance Hamiltonian. 

The n : m resonant Hamiltonian corresponding to the 
classical limit of the quantum Hamiltonian is: 

H(l, 6) = H (I) + 2tI™ ,2 I 2 1 ' 2 cos(m0i - n9 2 ) (19) 

The classical limit Hamiltonian can be reduced to an ef- 
fective one dimensional one by using the following gener- 
ating function 



f 2 = i ( 0t — e 2 ) + p c e 2 

171 



(20) 



The reduced Hamiltonian is obtained as 

H(I,<P;P C ) =H (I;P C ) +2rV(P,P c ) cos ™0 (21) 
with 

i/2 



V(I;P c )=I m ' 2 (P c --I 



(22) 
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The action P c is a constant of the motion since the 
Hamiltonian is ignorable in the angle conjugate to P c . 
P c = n/mli + I2 is the classical analog of the polyad 
number. 

The fixed points in the (/, <f) plane correspond to pe- 
riodic orbits in the full phase space. Using I = = <j) 
we obtain the fixed points as <fi = 0, ±7r/m. The solution 
I = is discarded since it is unphysical and the other 
root of I is of no interest here since that corresponds to 
all the excitation in the first mode. Using the value of 
the angle (f> at the fixed points the corresponding actions 
are obtained by a solution to the algebraic equation: 



dHo 
dl 



±2t 9 (I;P c ) = 



where 



g(i; Pc) 



dV(P,P c ) 
dl 



(23) 



(24) 



This nonlinear equation is solved numerically to obtain 
the periodic orbits of the system. 

In this appendix, however, attention is focused on the 
analysis of the function 



dH 



2V(I; P c ) cos mq 



(25) 



which occurs in the semiclassical approximation Eq. (jlj) 
for the quantum level velocities. Interest in this function 
arises from exploring the maximum and minimum values 
that this function can take at the fixed points in the (/, 4>) 



space. Inspecting the / and 4> derivative of v(I, <j)\ P c ) 



-v(i,<t>;Pc 



-2mV(I; P c ) shimi; 



01 



v(I,<P;P c ) = 2g(I;P c ) cos m<\> 



(26a) 
(26b) 



it is clear that the critical points in the <j> variable 
are identical with those of the fixed points. The 
critical points arising from the I variable are / = 
0,(m/n)P c ,m 2 P c /n(m + n). Note that demanding 
g(P, P c ) — is tantamount to being on the center of the 
resonance line and hence the stable periodic orbit. The 
first two critical values give rise to zero velocities and cor- 
respond to the "edge" states in a given polyad multiplet. 
The third critical point gives rise to the velocities: 



,(I (3 \0;P c ) = -v{f 3 \n/m;P c ) 

-1 1/2 



= 2 



(1 + M) 



n+rn 



(27a) 

p( m +n)/2 ( 27b ) 



where, /1 = n/m. From the above one immediately sees 
that the classical object takes the same absolute value at 
the stable and unstable fixed points. Moreover, the clas- 
sical analog exhibits the identical scaling with the polyad 
P c as seen in exact quantum studies done previously. In 
fact it is seen in the quantum studies that the resonant 
eigenstate \P;v) with v — i.e., state localized on the 
stable fixed point shows very good agreement with the 
above classical estimate. 
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